Systemic inflammation, the peripheral blood transcriptome and primary melanoma.

Peripheral blood transcriptomes from 383 newly-diagnosed melanoma patients were subjected to differential gene expression analysis. The hypotheses were that impaired systemic immunity is associated with poorer prognosis (thicker tumors and fewer tumor infiltrating lymphocytes (TILs)) and evidence of systemic inflammation (high-sensitivity C-reactive protein (hsCRP) and fibrinogen levels). Higher fibrinogen levels were associated with thicker primaries. In single gene analysis hsCRP levels were significantly associated with higher blood CD274 expression, (coding for PD-L1), but each was independently prognostic, with hsCRP associated with increased mortality, and higher CD274 protective, independent of age. Pathway analysis identified downregulation of immune cell signalling pathways in the blood of people with thicker tumors and notable upregulation of STAT1 in people with brisk TILs. Transcriptomic data provided evidence for increased NFB signalling with higher inflammatory markers but with reduction in expression of HLA class II molecules and higher CD274 suggesting that aberrant systemic inflammation is a significant mediator of reduced immune function in melanoma. In summary, transcriptomic data revealed evidence of reduced immune function in patients with thicker tumors and fewer TILs, at diagnosis. Inflammatory markers were associated with thicker primaries and independently with death from melanoma suggesting that systemic inflammation contributes to that reduced immune function.

Immune subtyping of melanoma whole slide images using multiple instance learning.

Determining early-stage prognostic markers and stratifying patients for effective treatment are two key challenges for improving outcomes for melanoma patients. Previous studies have used tumour transcriptome data to stratify patients into immune subgroups, which were associated with differential melanoma specific survival and potential predictive biomarkers. However, acquiring transcriptome data is a time-consuming and costly process. Moreover, it is not routinely used in the current clinical workflow. Here, we attempt to overcome this by developing deep learning models to classify gigapixel haematoxylin and eosin (H&E) stained pathology slides, which are well established in clinical workflows, into these immune subgroups. We systematically assess six different multiple instance learning (MIL) frameworks, using five different image resolutions and three different feature extraction methods. We show that pathology-specific self-supervised models using 10x resolution patches generate superior representations for the classification of immune subtypes. In addition, in a primary melanoma dataset, we achieve a mean area under the receiver operating characteristic curve (AUC) of 0.80 for classifying histopathology images into 'high' or 'low immune' subgroups and a mean AUC of 0.82 in an independent TCGA melanoma dataset. Furthermore, we show that these models are able to stratify patients into 'high' and 'low immune' subgroups with significantly different melanoma specific survival outcomes (log rank test, P< 0.005). We anticipate that MIL methods will allow us to find new biomarkers of high importance, act as a tool for clinicians to infer the immune landscape of tumours and stratify patients, without needing to carry out additional expensive genetic tests.

GABA Regulates Electrical Activity and Tumor Initiation in Melanoma.

Oncogenes can initiate tumors only in certain cellular contexts, which is referred to as oncogenic competence. In melanoma, whether cells in the microenvironment can endow such competence remains unclear. Using a combination of zebrafish transgenesis coupled with human tissues, we demonstrate that GABAergic signaling between keratinocytes and melanocytes promotes melanoma initiation by BRAFV600E. GABA is synthesized in melanoma cells, which then acts on GABA-A receptors in keratinocytes. Electron microscopy demonstrates specialized cell-cell junctions between keratinocytes and melanoma cells, and multielectrode array analysis shows that GABA acts to inhibit electrical activity in melanoma/keratinocyte cocultures. Genetic and pharmacologic perturbation of GABA synthesis abrogates melanoma initiation in vivo. These data suggest that GABAergic signaling across the skin microenvironment regulates the ability of oncogenes to initiate melanoma. SIGNIFICANCE: This study shows evidence of GABA-mediated regulation of electrical activity between melanoma cells and keratinocytes, providing a new mechanism by which the microenvironment promotes tumor initiation. This provides insights into the role of the skin microenvironment in early melanomas while identifying GABA as a potential therapeutic target in melanoma. See related commentary by Ceol, p. 2128. This article is featured in Selected Articles from This Issue, p. 2109.

Identifying the Transcriptional Drivers of Metastasis Embedded within Localized Melanoma.

In melanoma, predicting which tumors will ultimately metastasize guides treatment decisions. Transcriptional signatures of primary tumors have been utilized to predict metastasis, but which among these are driver or passenger events remains unclear. We used data from the adjuvant AVAST-M trial to identify a predictive gene signature in localized tumors that ultimately metastasized. Using a zebrafish model of primary melanoma, we interrogated the top genes from the AVAST-M signature in vivo. This identified GRAMD1B, a cholesterol transfer protein, as a bona fide metastasis suppressor, with a majority of knockout animals rapidly developing metastasis. Mechanistically, excess free cholesterol or its metabolite 27-hydroxycholesterol promotes invasiveness via activation of an AP-1 program, which is associated with increased metastasis in humans. Our data demonstrate that the transcriptional seeds of metastasis are embedded within localized tumors, suggesting that early targeting of these programs can be used to prevent metastatic relapse. SIGNIFICANCE: We analyzed human melanoma transcriptomics data to identify a gene signature predictive of metastasis. To rapidly test clinical signatures, we built a genetic metastasis platform in adult zebrafish and identified GRAMD1B as a suppressor of melanoma metastasis. GRAMD1B-associated cholesterol overload activates an AP-1 program to promote melanoma invasion. This article is highlighted in the In This Issue feature, p. 1.

ROR2 promotes epithelial-mesenchymal transition by hyperactivating ERK in melanoma.

Receptor tyrosine kinase-like orphan receptor 2 (ROR2) is a protein with important functions during embryogenesis that is dysregulated in human cancer. An intriguing feature of this receptor is that it plays opposite roles in different tumor types either promoting or inhibiting tumor progression. Understanding the complex role of this receptor requires a more profound exploration of both the altered biological and molecular mechanisms. Here, we describe that ROR2 promotes Epithelial-Mesenchymal Transition (EMT) by inducing cadherin switch and the upregulation of the transcription factors ZEB1, Twist, Slug, Snail, and HIF1A, together with a mesenchymal phenotype and increased migration. We show that ROR2 activates both p38 and ERK mitogen-activated protein kinase pathways independently of Wnt5a. Further, we demonstrated that the upregulation of EMT-related proteins depends on the hyperactivation of the ERK pathway far above the typical high constitutive activity observed in melanoma. In addition, ROR2 also promoted ERK phosphorylation, EMT, invasion, and necrosis in xenotransplanted mice. ROR2 also associates with EMT in tumor samples from melanoma patients where analysis of large cohorts revealed that increased ROR2 levels are linked to EMT signatures. This important role of ROR2 translates into melanoma patient' s prognosis since elevated ROR2 levels reduced overall survival and distant metastasis-free survival of patients with lymph node metastasis. In sum, these results demonstrate that ROR2 contributes to melanoma progression by inducing EMT and necrosis and can be an attractive therapeutic target for melanoma.

Association Study between Polymorphisms in DNA Methylation-Related Genes and Testicular Germ Cell Tumor Risk.

BACKGROUND: Testicular germ cell tumors (TGCT), histologically classified as seminomas and nonseminomas, are believed to arise from primordial gonocytes, with the maturation process blocked when they are subjected to DNA methylation reprogramming. SNPs in DNA methylation machinery and folate-dependent one-carbon metabolism genes have been postulated to influence the proper establishment of DNA methylation. METHODS: In this pathway-focused investigation, we evaluated the association between 273 selected tag SNPs from 28 DNA methylation-related genes and TGCT risk. We carried out association analysis at individual SNP and gene-based level using summary statistics from the Genome Wide Association Study meta-analysis recently conducted by the international Testicular Cancer Consortium on 10,156 TGCT cases and 179,683 controls. RESULTS: In individual SNP analyses, seven SNPs, four mapping within MTHFR, were associated with TGCT risk after correction for multiple testing (q ≤ 0.05). Queries of public databases showed that three of these SNPs were associated with MTHFR changes in enzymatic activity (rs1801133) or expression level in testis tissue (rs12121543, rs1476413). Gene-based analyses revealed MTHFR (q = 8.4 × 10-4), methyl-CpG-binding protein 2 (MECP2; q = 2 × 10-3), and ZBTB4 (q = 0.03) as the top TGCT-associated genes. Stratifying by tumor histology, four MTHFR SNPs were associated with seminoma. In gene-based analysis MTHFR was associated with risk of seminoma (q = 2.8 × 10-4), but not with nonseminomatous tumors (q = 0.22). CONCLUSIONS: Genetic variants within MTHFR, potentially having an impact on the DNA methylation pattern, are associated with TGCT risk. IMPACT: This finding suggests that TGCT pathogenesis could be associated with the folate cycle status, and this relation could be partly due to hereditary factors.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci.

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

Genes regulated by DNA methylation are involved in distinct phenotypes during melanoma progression and are prognostic factors for patients.

In addition to mutations, epigenetic alterations are important contributors to malignant transformation and tumor progression. The aim of this work was to identify epigenetic events in which promoter or gene body DNA methylation induces gene expression changes that drive melanocyte malignant transformation and metastasis. We previously developed a linear mouse model of melanoma progression consisting of spontaneously immortalized melanocytes, premalignant melanocytes, a nonmetastatic tumorigenic, and a metastatic cell line. Here, through the integrative analysis of methylome and transcriptome data, we identified the relationship between promoter and/or gene body DNA methylation alterations and gene expression in early, intermediate, and late stages of melanoma progression. We identified adenylate cyclase type 3 (Adcy3) and inositol polyphosphate 4-phosphatase type II (Inpp4b), which affect tumor growth and metastatic potential, respectively. Importantly, the gene expression and DNA methylation profiles found in this murine model of melanoma progression were correlated with available clinical data from large population-based primary melanoma cohorts, revealing potential prognostic markers.

Ulcerated melanoma: Systems biology evidence of inflammatory imbalance towards pro-tumourigenicity.

Microscopic ulceration is an independent predictor of melanoma death. Here, we used systems biology to query the role of host and tumour-specific processes in defining the phenotype. Albumin level as a measure of systemic inflammation was predictive of fewer tumour-infiltrating lymphocytes and poorer survival in the Leeds Melanoma Cohort. Ulcerated melanomas were thicker and more mitotically active (with corresponding transcriptomic upregulated cell cycle pathways). Sequencing identified tumoural p53 and APC mutations, and TUBB2B amplification as associated with the phenotype. Ulcerated tumours had perturbed expression of cytokine genes, consistent with protumourigenic inflammation and histological and transcriptomic evidence for reduced adaptive immune cell infiltration. Pathway/network analysis of multiomic data using neural networks highlighted a role for the ß-catenin pathway in the ulceration, linking genomic changes in the tumour to immunosuppression and cell proliferation. In summary, the data suggest that ulceration is in part associated with genomic changes but that host factors also predict melanoma death with evidence of reduced immune responses to the tumour.

ROR2 has a protective role in melanoma by inhibiting Akt activity, cell-cycle progression, and proliferation.

BACKGROUND: Receptor tyrosine kinase-like orphan receptor 2 (ROR2) is a Wnt5a receptor aberrantly expressed in cancer that was shown to either suppress or promote carcinogenesis in different tumor types. Our goal was to study the role of ROR2 in melanoma. METHODS: Gain and loss-of-function strategies were applied to study the biological function of ROR2 in melanoma. Proliferation assays, flow cytometry, and western blotting were used to evaluate cell proliferation and changes in expression levels of cell-cycle and proliferation markers. The role of ROR2 in tumor growth was assessed in xenotransplantation experiments followed by immunohistochemistry analysis of the tumors. The role of ROR2 in melanoma patients was assessed by analysis of clinical data from the Leeds Melanoma Cohort. RESULTS: Unlike previous findings describing ROR2 as an oncogene in melanoma, we describe that ROR2 prevents tumor growth by inhibiting cell-cycle progression and the proliferation of melanoma cells. The effect of ROR2 is mediated by inhibition of Akt phosphorylation and activity which, in turn, regulates the expression, phosphorylation, and localization of major cell-cycle regulators including cyclins (A, B, D, and E), CDK1, CDK4, RB, p21, and p27. Xenotransplantation experiments demonstrated that ROR2 also reduces proliferation in vivo, resulting in inhibition of tumor growth. In agreement with these findings, a higher ROR2 level favors thin and non-ulcerated primary melanomas with reduced mitotic rate and better prognosis. CONCLUSION: We conclude that the expression of ROR2 slows down the growth of primary tumors and contributes to prolonging melanoma survival. Our results demonstrate that ROR2 has a far more complex role than originally described.

Cooperation between melanoma cell states promotes metastasis through heterotypic cluster formation.

Melanomas can have multiple coexisting cell states, including proliferative (PRO) versus invasive (INV) subpopulations that represent a "go or grow" trade-off; however, how these populations interact is poorly understood. Using a combination of zebrafish modeling and analysis of patient samples, we show that INV and PRO cells form spatially structured heterotypic clusters and cooperate in the seeding of metastasis, maintaining cell state heterogeneity. INV cells adhere tightly to each other and form clusters with a rim of PRO cells. Intravital imaging demonstrated cooperation in which INV cells facilitate dissemination of less metastatic PRO cells. We identified the TFAP2 neural crest transcription factor as a master regulator of clustering and PRO/INV states. Isolation of clusters from patients with metastatic melanoma revealed a subset with heterotypic PRO-INV clusters. Our data suggest a framework for the co-existence of these two divergent cell populations, in which heterotypic clusters promote metastasis via cell-cell cooperation.

Identification of 22 susceptibility loci associated with testicular germ cell tumors.

Testicular germ cell tumors (TGCT) are the most common tumor in young white men and have a high heritability. In this study, the international Testicular Cancer Consortium assemble 10,156 and 179,683 men with and without TGCT, respectively, for a genome-wide association study. This meta-analysis identifies 22 TGCT susceptibility loci, bringing the total to 78, which account for 44% of disease heritability. Men with a polygenic risk score (PRS) in the 95th percentile have a 6.8-fold increased risk of TGCT compared to men with median scores. Among men with independent TGCT risk factors such as cryptorchidism, the PRS may guide screening decisions with the goal of reducing treatment-related complications causing long-term morbidity in survivors. These findings emphasize the interconnected nature of two known pathways that promote TGCT susceptibility: male germ cell development within its somatic niche and regulation of chromosomal division and structure, and implicate an additional biological pathway, mRNA translation.

Inherited duplications of PPP2R3B predispose to nevi and melanoma via a C21orf91-driven proliferative phenotype.

PURPOSE: Much of the heredity of melanoma remains unexplained. We sought predisposing germline copy-number variants using a rare disease approach. METHODS: Whole-genome copy-number findings in patients with melanoma predisposition syndrome congenital melanocytic nevus were extrapolated to a sporadic melanoma cohort. Functional effects of duplications in PPP2R3B were investigated using immunohistochemistry, transcriptomics, and stable inducible cellular models, themselves characterized using RNAseq, quantitative real-time polymerase chain reaction (qRT-PCR), reverse phase protein arrays, immunoblotting, RNA interference, immunocytochemistry, proliferation, and migration assays. RESULTS: We identify here a previously unreported genetic susceptibility to melanoma and melanocytic nevi, familial duplications of gene PPP2R3B. This encodes PR70, a regulatory unit of critical phosphatase PP2A. Duplications increase expression of PR70 in human nevus, and increased expression in melanoma tissue correlates with survival via a nonimmunological mechanism. PPP2R3B overexpression induces pigment cell switching toward proliferation and away from migration. Importantly, this is independent of the known microphthalmia-associated transcription factor (MITF)-controlled switch, instead driven by C21orf91. Finally, C21orf91 is demonstrated to be downstream of MITF as well as PR70. CONCLUSION: This work confirms the power of a rare disease approach, identifying a previously unreported copy-number change predisposing to melanocytic neoplasia, and discovers C21orf91 as a potentially targetable hub in the control of phenotype switching.

Transcriptional signatures underlying dynamic phenotypic switching and novel disease biomarkers in a linear cellular model of melanoma progression.

Despite advances in therapeutics, the progression of melanoma to metastasis still confers a poor outcome to patients. Nevertheless, there is a scarcity of biological models to understand cellular and molecular changes taking place along disease progression. Here, we characterized the transcriptome profiles of a multi-stage murine model of melanoma progression comprising a nontumorigenic melanocyte lineage (melan-a), premalignant melanocytes (4C), nonmetastatic (4C11-) and metastasis-prone (4C11+) melanoma cells. Clustering analyses have grouped the 4 cell lines according to their differentiated (melan-a and 4C11+) or undifferentiated/"mesenchymal-like" (4C and 4C11-) morphologies, suggesting dynamic gene expression patterns associated with the transition between these phenotypes. The cell plasticity observed in the murine melanoma progression model was corroborated by molecular markers described during stepwise human melanoma differentiation, as the differentiated cell lines in our model exhibit upregulation of transitory and melanocytic markers, whereas "mesenchymal-like" cells show increased expression of undifferentiated and neural crest-like markers. Sets of differentially expressed genes (DEGs) were detected at each transition step of tumor progression, and transcriptional signatures related to malignancy, metastasis and epithelial-to-mesenchymal transition were identified. Finally, DEGs were mapped to their human orthologs and evaluated in uni- and multivariate survival analyses using gene expression and clinical data of 703 drug-naïve primary melanoma patients, revealing several independent candidate prognostic markers. Altogether, these results provide novel insights into the molecular mechanisms underlying the phenotypic switch taking place during melanoma progression, reveal potential drug targets and prognostic biomarkers, and corroborate the translational relevance of this unique sequential model of melanoma progression.

MX2 mediates establishment of interferon response profile, regulates XAF1, and can sensitize melanoma cells to targeted therapy.

MX2 is an interferon inducible gene that is mostly known for its antiviral activity. We have previously demonstrated that MX2 is also associated with the tumorigenesis process in melanoma. However, it remains unknown which molecular mechanisms are regulated by MX2 in response to interferon signaling in this disease. Here, we report that MX2 is necessary for the establishment of an interferon-induced transcriptional profile partially through regulation of STAT1 phosphorylation and other interferon-related downstream factors, including proapoptotic tumor suppressor XAF1. MX2 and XAF1 expression tightly correlate in both cultured melanoma cell lines and in patient-derived primary and metastatic tumors, where they also are significantly related with survival. MX2 mediates IFN growth-inhibitory signals in both XAF1 dependent and independent ways and in a cell type and context-dependent manner. Higher MX2 expression renders melanoma cells more sensitive to targeted therapy drugs such as vemurafenib and trametinib; however, this effect is XAF1 independent. In summary, we uncovered a new mechanism in the complex regulation of interferon signaling in melanoma that can influence both survival and response to therapy.

Tumour gene expression signature in primary melanoma predicts long-term outcomes.

Adjuvant systemic therapies are now routinely used following resection of stage III melanoma, however accurate prognostic information is needed to better stratify patients. We use differential expression analyses of primary tumours from 204 RNA-sequenced melanomas within a large adjuvant trial, identifying a 121 metastasis-associated gene signature. This signature strongly associated with progression-free (HR = 1.63, p = 5.24 × 10-5) and overall survival (HR = 1.61, p = 1.67 × 10-4), was validated in 175 regional lymph nodes metastasis as well as two externally ascertained datasets. The machine learning classification models trained using the signature genes performed significantly better in predicting metastases than models trained with clinical covariates (pAUROC = 7.03 × 10-4), or published prognostic signatures (pAUROC < 0.05). The signature score negatively correlated with measures of immune cell infiltration (ρ = -0.75, p < 2.2 × 10-16), with a higher score representing reduced lymphocyte infiltration and a higher 5-year risk of death in stage II melanoma. Our expression signature identifies melanoma patients at higher risk of metastases and warrants further evaluation in adjuvant clinical trials.

Application of Circulating Cell-Free Tumor DNA Profiles for Therapeutic Monitoring and Outcome Prediction in Genetically Heterogeneous Metastatic Melanoma.

PURPOSE: Circulating cell-free tumor DNA (ctDNA) reflects the heterogeneous spectrum of tumor-specific mutations, especially in systemic disease. We validated plasma-based assays that allow the dynamic quantitative detection of ctDNA as a prognostic biomarker for tumor load and prediction of therapy response in melanoma. MATERIALS AND METHODS: We analyzed plasma-derived ctDNA from a large training cohort (n = 96) of patients with advanced-stage melanoma, with assays for the BRAF V600E and NRAS Q61 driver mutations as well as TERT C250T and TERT C228T promoter mutations. An independent patient cohort (n = 35) was used to validate the utility of ctDNA monitoring under mitogen-activated protein kinase-targeted or immune checkpoint therapies. RESULTS: Elevated plasma ctDNA level at baseline was an independent prognostic factor of disease progression when compared with serum S100 and lactate dehydrogenase levels in multivariable analyses (hazard ratio [HR], 7.43; 95% CI, 1.01 to 55.19; P = .05). The change in ctDNA levels during therapy correlated with treatment response, where increasing ctDNA was predictive for shorter progression-free survival (eg, for BRAF V600E ctDNA, HR, 3.70; 95% CI, 1.86 to 7.34; P < .001). Increasing ctDNA levels predicted disease progression significantly earlier than did routine radiologic scans (P < .05), with a mean lead time of 3.5 months. NRAS-mutant ctDNA was detected in a significant proportion of patients with BRAF-mutant tumors under therapy, but unexpectedly also at baseline. In vitro sensitivity studies suggested that this represents higher-than-expected intratumoral heterogeneity. The detection of NRAS Q61 ctDNA in baseline samples of patients with BRAF V600E mutation who were treated with mitogen-activated protein kinase inhibitors significantly correlated with shorter progression-free survival (HR, 3.18; 95% CI, 1.31 to 7.68; P = .03) and shorter overall survival (HR, 4.08; 95% CI, 1.57 to 10.58; P = .01). CONCLUSION: Our results show the potential role of ctDNA measurement as a sensitive monitoring and prediction tool for the early assessment of disease progression and therapeutic response in patients with metastatic melanoma.

Gene co-expression and histone modification signatures are associated with melanoma progression, epithelial-to-mesenchymal transition, and metastasis.

BACKGROUND: We have previously developed a murine cellular system that models the transformation from melanocytes to metastatic melanoma cells. This model was established by cycles of anchorage impediment of melanocytes and consists of four cell lines: differentiated melanocytes (melan-a), pre-malignant melanocytes (4C), malignant (4C11-), and metastasis-prone (4C11+) melanoma cells. Here, we searched for transcriptional and epigenetic signatures associated with melanoma progression and metastasis by performing a gene co-expression analysis of transcriptome data and a mass-spectrometry-based profiling of histone modifications in this model. RESULTS: Eighteen modules of co-expressed genes were identified, and some of them were associated with melanoma progression, epithelial-to-mesenchymal transition (EMT), and metastasis. The genes in these modules participate in biological processes like focal adhesion, cell migration, extracellular matrix organization, endocytosis, cell cycle, DNA repair, protein ubiquitination, and autophagy. Modules and hub signatures related to EMT and metastasis (turquoise, green yellow, and yellow) were significantly enriched in genes associated to patient survival in two independent melanoma cohorts (TCGA and Leeds), suggesting they could be sources of novel prognostic biomarkers. Clusters of histone modifications were also linked to melanoma progression, EMT, and metastasis. Reduced levels of H4K5ac and H4K8ac marks were seen in the pre-malignant and tumorigenic cell lines, whereas the methylation patterns of H3K4, H3K56, and H4K20 were related to EMT. Moreover, the metastatic 4C11+ cell line showed higher H3K9me2 and H3K36me3 methylation, lower H3K18me1, H3K23me1, H3K79me2, and H3K36me2 marks and, in agreement, downregulation of the H3K36me2 methyltransferase Nsd1. CONCLUSIONS: We uncovered transcriptional and histone modification signatures that may be molecular events driving melanoma progression and metastasis, which can aid in the identification of novel prognostic genes and drug targets for treating the disease.

Study of the Female Sex Survival Advantage in Melanoma-A Focus on X-Linked Epigenetic Regulators and Immune Responses in Two Cohorts.

BACKGROUND: Survival from melanoma is strongly related to patient sex, with females having a survival rate almost twice that of males. Many explanations have been proposed but have not withstood critical scrutiny. Prior analysis of different cancers with a sex bias has identified six X-linked genes that escape X chromosome inactivation in females and are, therefore, potentially involved in sex differences in survival. Four of the genes are well-known epigenetic regulators that are known to influence the expression of hundreds of other genes and signaling pathways in cancer. METHODS: Survival and interaction analysis were performed on the skin cutaneous melanoma (SKCM) cohort in The Cancer Genome Atlas (TCGA), comparing high vs. low expression of KDM6A, ATRX, KDM5C, and DDX3X. The Leeds melanoma cohort (LMC) on 678 patients with primary melanoma was used as a validation cohort. RESULTS: Analysis of TCGA data revealed that two of these genes-KDM6A and ATRX-were associated with improved survival from melanoma. Tumoral KDM6A was expressed at higher levels in females and was associated with inferred lymphoid infiltration into melanoma. Gene set analysis of high KDM6A showed strong associations with immune responses and downregulation of genes associated with Myc and other oncogenic pathways. The LMC analysis confirmed the prognostic significance of KDM6A and its interaction with EZH2 but also revealed the expression of KDM5C and DDX3X to be prognostically significant. The analysis also confirmed a partial correlation of KDM6A with immune tumor infiltrates. CONCLUSION: When considered together, the results from these two series are consistent with the involvement of X-linked epigenetic regulators in the improved survival of females from melanoma. The identification of gene signatures associated with their expression presents insights into the development of new treatment initiatives but provides a basis for exploration in future studies.